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Abstract 



We study the stability of the Higgs potential in the framework of the effective 
Lagrangian beyond the MSSM. While the leading nonrenormalizable operators can 
shift the Higgs boson mass above the experimental bound, they also tend to render 
the scalar potential unbounded from below. The destabilization is correlated with the 
Higgs mass increase, so that if quantum corrections are small the problem is severe. We 
show that a supersymmetric sub-leading correction stabilizes the potential within the 
domain of validity of the effective theory. Constraints on MSSM parameters as well as 
on higher dimensional operators are derived, ensuring that our vacuum has a lifetime 
longer than the present age of the universe. In addition we show that when effective 
operators are responsible for evading the LEP bound, stability constraints imply an 
upper bound on the scale of new physics in the few TeV range. 



1 Introduction 



The Higgs tree level quartic couplings in the Minimal Supersymmetric Standard Model 
(MSSM) are completely dictated by gauge interactions. This fact stands at the heart of 
theoretical and phenomenological difficulties, perhaps foremost among them the tree level 
prediction that the lightest Higgs boson be lighter than the Z, a prediction ruled out by 
collider experiments [HE]- In order to evade the experimental bound on the Higgs mass large 
quantum corrections are required in the MSSM, implying a substantial hierarchy between the 
electroweak (EW) scale and the scale of supersymmetry (SUSY) breaking. In particular, at 
least one of the top superpartners is required to be much heavier than the top such that some 
amount of fine-tuning is needed [31 HI El E] • The fine-tuning becomes even more pronounced 
if the model is expected to provide explanations for present cosmological data. While the 
MSSM possesses all of the ingredients needed in order to account for both the dark matter 
and the baryon asymmetry of the universe, analyses reveal that difficulties associated with 
the Higgs mass bound are rooted in the cosmological arcade, too (see, for example [HIEIE])- 

The restricted structure of the Higgs sector also makes it susceptible to small corrections 
from new physics Beyond the MSSM (BMSSM). If the scale associated with the BMSSM 
physics lies well above MSSM particle masses, an effective field theory approach becomes 
useful. The effective theory framework of the BMSSM was previously studied, e.g., in [TOlfTTl 
H2]. The authors of [TT] showed that the effective expansion takes on a rather simple form. In 
fact, under mild assumptions the leading nonrenormalizable corrections to the Higgs sector 
are captured by only two operators, one supersymmetric and the other associated with hard 
SUSY breaking (see also [IS])- It was further demonstrated that these operators may lead to a 
sizable shift of the Higgs mass at the classical level. Besides potentially solving the difficulties 
mentioned above [El [El [16], this result is of considerable experimental importance as it 
opens up a zone of SUSY phenomenology in which both stops may be very light, possibly 
just around the corner for collider experiments. 

However, examining the vacuum structure of the effective theory reveals that when the 
leading nonrenormalizable operators are taken to account for a significant shift for the lightest 
Higgs mass, they also destabilize the quartic couplings. Naively, this might be considered a 
severe setback to the picture drawn above. In the presence of a negative quartic coupling, 
stability relies on the higher order terms of the theory, which could a priori complicate 
the analysis considerably as well as introduce UV sensitivity. Nevertheless, we find that 
the effective theory exhibits a remarkable property: a single higher order operator arising 
at dimension six automatically cures the runaway initiated by the leading dimension five 
terms. Under mild assumptions, this operator is the only one relevant to the question of 
stability, and, furthermore, is directly correlated with the dimension five set. We find that 
the entire study of potential stability can be conducted and resolved within the range of 
validity of the effective approach. As a result, vacuum stability can be ensured by imposing 
simple relations between MSSM parameters and dimension five BMSSM operators, without 
the need to explicitly invoke the potentially far more complicated dimension six structure. 
In addition, we show that an upper bound on the scale of new physics may arise in the 
BMSSM, as well as derive the viable range for the lightest Higgs mass. 

The interplay between the negative quartic and the higher dimensional operator often 
leads to the formation of a remote vacuum, in the presence of which the EW vacuum is only 
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metastable. The existence of the remote vacuum was noticed by the authors of [T7j, where 
it was suggested that EW symmetry breaking may in fact lead our universe to reside in this 
new configuration dominated by nonrenormalizable operators. An intriguing feature of this 
scenario is that in this case, EW breaking may occur even in the SUSY limit. However, 
the structure of the effective theory implies that the EW scale be given in this case by the 
geometric mean of the new physics scale M and the SUSY fj, parametei0, v ~ yj r jlM. Null 
searches for charginos, for instance, highly restrict this possibility. In this paper we adopt 
a more conservative approach. We discuss the possibility that while the leading effective 
operators act to destabilize the potential, the usual EW breaking vacuum is either stable or 
of a life time longer than the present age of the universe. 

The outline of this paper is as follows. In Section [2] we review the BMSSM Higgs sector, 
emphasizing the correlation between the lightest Higgs mass and the appearance of directions 
in field space where the potential is unbounded from below at leading order. We then 
study the next set of higher dimensional operators, showing that they can stabilize the 
potential up to the cutoff scale of the effective theory. In Section [3] we present analytical 
and numerical prescriptions ensuring the (meta) stability of the EW vacuum. Section H] 
discusses the phenomenological implications of the stability constraints. Our conclusions are 
gathered in Section El Appendix |A] addresses the issue of charge breaking and CP violating 
field configurations, Appendix [B] elaborates on the quantum tunneling computation, and 
Appendix ICl presents constraints on the heavy cutoff scale arising from electroweak precision 
measurements. 

2 Higgs sector in the MSSM and Beyond 

In the bulk of this paper we analyze the vacuum stability of the BMSSM theory, namely 
the MSSM Higgs sector augmented by nonrenormalizable operators. It is useful to adopt 
the nomenclature of [UJ and to classify the nonrenormalizable operators according to their 
scaling dimension, which is the total mass dimension of the fields contained in an operator, 
and their effective dimension, counting the powers of 1/M which suppress it. For instance, 
an operator such as (fi/M)h u hd\h u \ 2 +h.c. is of scaling dimension four but effective dimension 
five. 

In Section 12.11 we briefly discuss the renormalizable MSSM Higgs sector, and highlight 
the main features which guide us in the study of nonrenormalizable corrections to it. In 
Section 12.21 we review BMSSM operators of effective dimension five. These operators alone 
suffice to lift the lightest Higgs mass above the LEP bound. However, as we show, these 
operators can destabilize the scalar potential at large field values which are still within the 
domain of validity of the effective theory. In that case, naively, the effective theory truncated 
at this order is not consistent; higher order operators must cure the instability. The next set 
of operators consists of effective dimension six, and we study it in Section 12.31 

lr The emergence of an intermediate scale of similar formal form was previously considered in [18] . where 
it was discussed in the context of much higher energy phenomena. 
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2.1 MSSM setup 

The scalar Higgs potential of the renormalizable MSSM can be written as followfl 
Vmssm = rn\ \h d \ 2 + ml \h u \ 2 + {m\ 2 h u h d + h.c.) 




where m\ = m 2 H + \p\ 2 , m\ = m 2 Hu + \p\ 2 , m\ 2 = Bp and g 2 z = g 2 +g' 2 with g, g' the SM gauge 
couplings. One can always choose a basis for the fields such that m\ 2 is real and positive, 
and we shall keep to such a basis consistently throughout the paper. We parameterize the 
expectation values as follows: 

(M T = (O,0 2 ), (h d ) T = {<f>i + ix,p)- (2) 

Gauge freedom is used to render = and 02 real and positive, while the real component 
0i is allowed to obtain negative valuesjfl In the BMSSM, as we shall see, more than one 
vacuum configuration may develop. We will analyze paths of minimum potential energy 
which connect these vacua. Along such paths, in principle, it may become energetically 
favorable for non-vanishing values of x or P t° turn on. This may occur even when all 
Lagrangian parameters are real, as we shall assume in this paper. In this case, we have 
found that ignoring the CP violating (CPV) and charge-breaking (CB) background fields x 
and p is typically well justified in the study of potential stability. Therefore, for clarity, we 
set p = x — in most of the paper. In places where deviations from this assumption become 
relevant we refer the reader to the discussion in Appendix [A] 

In the framework we study, SU(2)l x U(1)y breaking into U(1)em occurs as usual by an 
interplay between the quadratic and quartic terms in the scalar potential. (For an alternative 
scenario, see |17j.) The resulting vacuum is parameterized by 

tan/5 = 2 /0i. (3) 

We call this vacuum the EW vacuum. 

Being exclusively dictated by gauge superfield D-terms, the quartic couplings in (TfJ) are 
proportional to the electroweak gauge couplings, making the Higgs sector sensitive to small 
corrections. In particular, quantum corrections arising from loops of MSSM fields with a 
large amount of SUSY breaking are usually conceived to account for a sizable shift in the 
lightest Higgs mass. In a similar manner, the quartic Higgs structure is also sensitive to 
nonrenormalizable operators, which need only compete with couplings of order g 2 in order 
to modify the spectrum. 

Moreover, the quartic terms of §T§ vanish along the D-flat directions, specified by | (h u ) \ = 
| (ha) |. At tree (and renormalizable) level, this gives rise to a constraint on the quadratic 

2 We write the superfield components as H = (h,tp,F), the MSSM /j, term as f d 2 8iiH u Hd, and our 
convention for SU(2) contraction is H u Hd = H+H^ — H^H®. 

3 Note that, in the renormalizable MSSM, the minimization conditions lead to p = \ = 0. In that case 
one may write the scalar potential using non- negative <f>\ and 02 . 
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terms, m\ + m\ > 2m 2 2 , which must hold for the potential to be bounded from below. 
When nonrenormalizable operators are considered, the relative flatness of the renormalizable 
potential along the D-flat directions makes it important to verify that the higher dimensional 
operators do not destabilize the vacuum. 



2.2 Operators of effective dimension five 

Effective dimension five operators composed purely of Higgs fields enter the Lagrangian via 
the superpotential. Including F-term SUSY breaking, there are two such operators [TT] : 

J d 2 6 {\ x {H u H d f + X 2 Z(H u H d ) 2 ) + h.c. (4) 

Here Z = ms 9 2 is a dimensionless chiral superfield spurion. One could also contemplate the 
existence of operators arising from D-term SUSY breaking, in which case additional effective 
dimension five operators arise. In this paper, however, we assume that the effect of D-term 
breaking is somehow suppressed or non-existent [TTl [T9] . 

We assume that the new physics generating the effective operators is approximately 
super symmetric, ms <C M. We are interested in the imprint which the effective dimension 
five operators have on both the spectrum and the stability of the scalar potential, as we now 
discuss. At effective dimension five, the correction to the scalar potential resulting from (Tjfl) 
reads: 

5V 5 = 2e x h u h d {\h u \ 2 + \h d \ 2 ) + e 2 (h u h d ) 2 + h.c, (5) 

where we have defined t\ = Ai/i*/M and e 2 = —X 2 ms/M. Expanding to order 0(e), the 
following shift is obtained for the light (CP-even) Higgs boson mass: 

2 , / 2e lr (m 2 A + m|) sin 2/3 + e 2r (m 2 A - m|) cos 2 2/3 \ 
S e m h = 2v e 2r -2e lr sm2/3 ========= == , (6) 

V V \ m A ~ m z) + 4m^m| sin 2 2(3 J 

where e& r denotes the real part of e^. As explained in [TT], only the real parts of enter the 
spectrum at leading order while the imaginary parts contribute to interactions and mixing. 
In this paper we make use of the fact that the spectrum is relatively insensitive to the 
imaginary parts of £1,2 and consider, for simplicity, only the case where there exists a basis 
for h u , h d in which €\ t2 are real and m\ 2 is real and positive. Henceforth we drop the r 
subscript and assume that ei i2 are real. 

The contribution of the non-supersymmetric term to the mass shift is suppressed com- 
pared to the supersymmetric one. For example, in the limit m\ ^> m 2 z we have 

ml 9 „ 4e 2 sin 2 2/3 16eisin2/3 / ' m 2 7 \ 

-4 = cos 2 2(3 + 2 H ^ — - + [-§■). 7 

m z 3z 3z \ m A/ 

It follows that IS^ml/ 5 ei m 2 \ ~ |e 2 /ei | sin 2/5 /4, and so e 2 > 4|ei| is needed in order for 
both terms to give a comparable mass shift. In Figure [JJ we report the size of both the 
SUSY preserving and breaking operators, as required for m^ = 115 GeV, illustrating that 
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Figure 1: Contour plots of the values of —e\ (left) and e<i (right) corresponding to a fixed 
Higgs mass rrih = 115 GeV. In each panel only the considered operator is nonzero. 



the former can easily lift the Higgs mass classically at this order, while existence of only the 
latter calls for large quantum corrections as in the MSSM. For moderate tan (3 a negative e\ 
of magnitude |ei| > 10~ 2 X (1 — r -2 ) tan/?, with r = m^/m^, can lift the Higgs mass above 
the LEP bound at tree level, without the need for quantum corrections. The large tan/3 
limit is more involved. For tan/3 > |l/ei|, the leading contribution of effective dimension 
five operators is suppressed to the level of the next order in e and, as a result, becomes 
comparable to that of effective dimension six terms. 

Apart from affecting the local properties of the vacuum {e.g. particle spectrum and 
interactions) the effective dimension five operators also influence the potential at large field 
values. Along the D-flat directions^ the potential is only bounded from below as long as 
4|ei| < £2- If this relation does not hold, the potential eventually becomes unstable and a 
runaway occurs along the direction </>i = sign (ej) 02- In the interesting case where e\ < 0, 
corresponding to a positive shift to the lightest Higgs mass, a saddle point emerges at field 
values 

\ m\ + mj + 2^ 

^ , = * 2 = V 8(| £l |-i £2 ) • (8) 

after which the potential decreases indefinitely. An important property of the runaway 
is that it occurs at field values </> 2 ~ m 2 /e ~ mM, where m refers to a combination of 
relevant quadratic mass parameters (presumably of electroweak scale) and M is the heavy 
scale. Therefore, the runaway develops at <p M, rendering the EW vacuum unacceptably 
short-lived. 

4 Note that there are now two such directions as in our convention (f>i can be negative. 
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Given the discussion above, one might naively conclude that the inequality 4|ei| < t2 
must be enforced as a physical constraint on the relative magnitudes of the supersymmetric 
vs. the SUSY breaking corrections [TJ]. In particular, in the supersymmetric limit wherein 
62 = 0, also ei = would be required. However, as we show below, the supposed runaway 
is an artifact of the truncation of the effective expansion at 0(M~ l ) [20] • This expansion 
cannot be trusted at large field values, even though it is consistent in the local vicinity of our 
vacuum [T7]. Indeed, effective dimension six operators generate positive scaling dimension 
six contributions of the form <j) e /M 2 . Such terms eventually win over the negative quartic 
terms at field values <fi 2 ~ mM, precisely the region of field space at which the instability 
starts to develop. Thus, effective dimension six operators are as important as the effective 
dimension five ones around the instability and may stabilize the potential before the UV 
threshold of the effective theory. Below we inspect the structure of effective dimension six 
operators in order to determine which ones are relevant to the stability problem. Compared 
to previous studies (2TJ [15] , we find that focusing on the physical problem at hand simplifies 
the analysis considerably. Delightfully, we are able to show that over a broad range of 
parameters only one such operator exists, and even this operator not strictly independent 
from the dimension five set. 



2.3 Operators of effective dimension six 

We begin by considering effective dimension six operators which involve SUSY breaking. At 
order 1/M 2 , there are no scaling dimension six operators which arise from SUSY breaking 
terms in the Lagrangian. The reason is that each term arising from SUSY breaking spurion 
must be accompanied by the appropriate power of the SUSY breaking mass scale. At scaling 
dimension six, this implies that the first (pure Higgs) operators are of order 1/M 3 . 

We do find effective operators associated with SUSY breaking spurions at scaling dimen- 
sion four. These are of the form (m/M) 2 h A , where we extend the definition of m to include 
the SUSY breaking scale, also of order the electroweak scale in this framework. Along 
the D-flat directions such contributions are suppressed by an additional power of m/M in 
comparison to the effective dimension five terms of 6V5. Away from the D-flat directions, 
MSSM D-terms guarantee stability provided that we impose 

~ 8 15 1 } 

as a parametric inequality, where e 2 = (m/M) 2 represents the magnitude of the dimensionless 
coefficient of such operators. Henceforth we restrict the discussion to the scenario in which 
Eq. ([HI) is satisfied. We shall see that this assumption simplifies the problem considerably. 
Moreover, it stands in accordance with requiring the effective theory expansion in powers 
of e to remain valid, regardless of the stability analysis. Referring to Figure [T] we find that, 
at least for moderate tan/?, Eq. @ does not pose any real limitation on the role of the ei 
correction in lifting the Higgs mass above the current experimental bound. Nevertheless, 
in general, the scaling dimension four, effective dimension six operators do not display any 
particular tan (3 dependence and so they can have a comparable influence on the spectrum in 
the large tan (3 regime, where the leading dimension five contributions are suppressed. Due 
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to the multiplicity of independent coefficients, the analysis of the spectrum is involved in 
this case and we do not pursue it further in this paper. 

From the discussion above we conclude that it is enough for our purpose to study effective 
dimension six operators in the supersymmetric limit. Squaring the 1/M piece of the F-term 
equations of motion, arising in the presence of the supersymmetric effective dimension five 
operator gives 

2 

\h u h d \ 2 (\h u \ 2 + \h d \ 2 ). (10) 

The contribution ( TTOT) is positive definite and non- vanishing along the D-flat directions. Thus 
it plays an important role in stabilizing the potential. Note that, for a given value of the 
H parameter, the coefficient of this contribution is correlated with the Higgs mass through 
Eq. ®. 

Superpotential operators which involve gauge superfields contribute to the scalar poten- 
tial at effective dimension six, and must contain D-term components. Kahler operators which 
arise at order 1/M 2 can only affect the scalar potential through either F-terms or D-terms, 
by gauge invariance. Considering F-terms, we immediately see that these cannot contribute 
to effective dimension six, scaling dimension six operators, since they are linear in h at 
leading order in 1/M. F-terms do contribute to effective dimension six, scaling dimension 
four operators. However, following Eq. (jHJ) and the related discussion, such contributions 
are parameterically suppressed and we need not pursue them further. Considering D-terms, 
and ignoring sub-leading scaling dimension four operators, we find scaling dimension six 
operators which are of the general form 

|UW-N 2 ) (ii) 

where g 2 stands for some bilinear combination of g and g' . 

We learn that, apart from the operator ffTOj) . all other scaling dimension six contributions 
to the scalar potential are gauge coupling suppressed and - most importantly - vanish along 
the D-flat directions of the MSSM. Therefore at large field values (but still <C M) and 
along the would-be runaway direction described in the previous section, the potential is in 
fact driven by the positive definite scaling dimension six SVq given in (TlOl) . Thus the same 
superpotential operator responsible for lifting the Higgs mass classically also ensures that 
the scalar potential does not exhibit a runaway. The phenomenology of this correlation we 
aim to address in the next section. 

3 Vacuum stability 

In what follows we analyze the stability of the scalar potential including the dimension 
six operator identified in the previous section, namely V = Vmssm + 8V5 + SVq. Simple 
analytical and numerical criteria which ensure the stability of a given potential configuration 
are formulated. These criteria can be put in terms of relations between the electroweak scale 
parameters and effective dimension five operators. 
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3.1 Analytical approximation 

It is useful to first analyze the potential along the MSSM D-flat directions. While soft 
terms, quantum corrections and the presence of the EW vacuum all play a role in shifting 
the potential features somewhat away from \{h u )\ = \ (hd)\, all of the insights are contained 
and, further more, it turns out to be a reasonable approximation to study the profile of the 
potential at these well defined directions in field space. We begin by performing this analysis 
at tree level, assuming vanishing \ and p values. Then, having obtained the principal results 
we extend the discussion to include all of the complications mentioned above. 

At tree level, the effective potential along the MSSM D-flat directions take the form: 



V D ~ flat (0) = \ {ml + m\ T 2m 2 2 ) 2 + 2 (| T e 1 



)0 4 + 









:,6 



(12) 



with ±0! = 02 = 4>/V%- If the quartic coupling in (TT2]) is negative, the potential may 
develop another vacuum away from the electroweak scale. Thus a simple way to guarantee 
stability is to impose €2 > 4|ei|. However, from Figure [TJwe learn that typical values for |ei|, 
consistent with the bound on the lightest Higgs mass are in the range |ei| ~ 0.05 — 0.1. The 
simple condition for an always-positive quartic is therefore in some tension with the need to 
assure that higher order terms are under control (see Eq. ©)• Hence we attend to the more 
interesting case where the quartic coupling is negative. It is useful to define the quantities 
m 2 = ml + m 2 =F 2m 2 2 and e = €2/4 =)= e\. Using these quantities we see that the non-zero 
extrema of the potential, corresponding to the remote vacuum and saddle point are located 
respectively at 

,2 2|/i| 2 /e 



3e Vei 




l±Wl-ST7Ti T • (13) 



Manipulating Eqs. (jT2"j) and (I13p we arrive at the following criterion, designed to ensure that 
the remote minimum along the D-flat direction is at most degenerate with, but never deeper 
than the potential at the origin of field space: 



m 2 I e 



w £ 2 UJ • (14) 

A remarkable feature of the condition CHI) is that it does not involve the heavy scale 
M. Rather, it involves a relation between relevant mass parameters of the MSSM and the 
relative size of nonrenormalizable operators. Indeed, from Eq. ( fl3l) it follows that both the 
remote vacuum and the saddle point scale similarly, as <fi ~ \JmM M. This behavior 
reflects the fact that the quartic couplings along the D-flat directions originate from effective 
dimension five operators. As a result, there exists a scale at which quadratic, quartic and 
dimension six contributions in the potential are of similar magnitude; this scale is precisely 
VmM. 

For concreteness we attend to the case €\ < 0, corresponding to a positive shift for the 
lightest Higgs mass. Rewriting Eq. (|14p in terms of physical quantities and to leading order 
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in e, we arrive at the following useful relation: 



1 + 



rrir 



16e /l + 2sin2/3 



m , 



g\ \ 1 + sin 2{3 



3ei 
2 e 



-i 



(15) 



Restricting to the D-flat direction allowed us to write a stability criterion of simple analytical 
form. However, the actual vacua and saddle point emerge somewhat away from the D-flat 
direction. In particular, the remote vacuum (and, similarly, saddle point) can sustain an 
angle ~ m/ymM ~ sfe from the flat direction. Hence imposing stability along |0i| = |02| 
does not forbid the actual remote vacuum from becoming a global minimum. It is therefore 
important to complement Eq. (fT5l) using numerical methods. In practice, as we show below, 
Eq. fll5p turns out to be a robust but slightly conservative stability criterion. 



3.2 Numerical approach 

Armed with intuition from the analytical analysis, we proceed to define stability criteria 
based on numerical procedures. We formulate two such criteria, then comment on the 
possibility of CPV or CB field configurations. 



Vacuum degeneracy First, we define a stability criterion by computing the potential 
and numerically verifying that the remote vacuum is at most degenerate with, but never 
deeper than the EW vacuum. We call this constraint vacuum degeneracy; it is robust but 
conservative. 



Quantum tunneling Second, we define a stability criterion by numerically computing 
the tunneling rate from the EW to the remote vacuum. The tunneling rate is given by 
T oc exp (— B), with B the bounce action. A metastable configuration is viable if the 
universe remains in the EW vacuum for longer than its age. Quantitatively, this translates 
into B > 400. We have used an approximate method to compute the bounce action, and 
errors of 0(1) are expected. Hence we discuss configurations with bounce action of B = 400 
and 10 3 . A detailed computation of the tunneling rate can be found in Appendix [Bi 



In the discussion above we have ignored the CP violating and charge breaking background 
fields, x an d P- However, these effects are accounted for when we numerically search for 
vacuum degeneracy and compute the bounce action. We find that neglecting the x an d 
p background fields is almost always justified in the analysis of potential stability in our 
framework, and does not affect the reliability of Eq. ( JT5il . The reason is that the stability 
criterion itself is typically sufficient to ban non-vanishing x an d p, for field values <fi < M. 
More details can be found in Appendix [A] 

In Figure [5] (three panels on the left) the stability criteria are projected on various sections 
of the parameter space, illustrating the above analysis. We focus on the combined study of 
electroweak scale parameters vs. the heavy scale M. The following statements can be made. 

The analytical stability criterion is slightly less conservative than, but closely follows the 
robust numerical criterion of vacuum degeneracy. Moreover, in the parameter space depicted 
in Figure [2] we find that the analytical criterion is typically more conservative than imposing 
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B > 10 3 , and always ensures that the bounce action will be above B ~ 400. Thus we can 
safely say that the analytical criterion is robust, at least for moderate values of tan/3 < 10 
where small ei < 0.1 is more than sufficient to lift the Higgs mass above the experimental 
bound. Since the gap between the various criteria is modest at all parameter values, it is 
evident that Eq. (fl5l) provides a detailed qualitative, as well as quantitative understanding 
of the parameter space of vacuum stability. 

Regarding the tunneling action, we find that it rises steeply above B ~ 100, rendering the 
B = 400 and 10 3 contours very close to each other. This occurs since in the relevant regions 
of parameter space the remote minimum is nearly degenerate with the EW vacuum, and so 
a small change of parameters towards vacuum degeneracy causes the action to diverge. 

4 Phenomenological implications 

Finally the stage is set to study the implications of stability constraints. We begin by 
discussing relations between electroweak scale parameters and the effective dimension five 
operators. We then study constraints on the heavy BMSSM scale M. Lastly we describe 
the viable shift to the lightest Higgs mass, and outline how the analysis is altered in the 
presence of radiative corrections. 

4.1 Constraints on BMSSM parameters 

The parameter space consistent with vacuum stability, depicted in Figure [5] for rrih = 115 
GeV, becomes smaller when the value of rrih is increased further above the experimental 
boundEI Hence Figure |2] can be interpreted as a translation of our current knowledge, rrih > 
114 GeV, into bounds on the viable ranges of other BMSSM parameters, subject to the 
assumption that quantum corrections are small. Furthermore, one of the merits of the 
analytical criterion Eq. ffT5"]) is that it allows to extract limiting relations between BMSSM 
parameters regardless of the precise value of the Higgs mass, as long as the expansion to 
order 0(e) is valid in the local vicinity of the EW vacuum. In order to gain further insight 
it is useful to look at the behavior of Eq. (fl5l) in the following two limits. 

Large rriA limit — m\ ^> m\: In this case Eq. ( fl5l) implies 



This condition is robust but conservative. In other words, while having < mA\ti/i\ 
ensures stability, somewhat larger values of \ii\/mA may still yield stable vacuum configura- 
tions, as illustrated in Figure [21 Indeed, comparing the two upper panels we see that while 
the vacuum degeneracy criterion requires > \fi\, the more loose tunneling action criterion 
allows rriA ~ \ft\. For example, taking tan/3 = 4, ei = — £i = 0.06, \fi\ = rriA = 3mz gives 

5 A quick way to understand this statement (which will be assessed in Section 1*0]) is by noting that, for 
given values of rriA, tan (3 and 62, the way to increase rrih is via making ei more negative. Vacuum instability 
is then driven stronger by the negative quartic ~ £i</> 4 , balanced at large fields by the dimension six term 
~ <j) 6 /M 2 . It follows that the upper bound on M derived from stability must decrease as rrih is increased. 
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Figure 2: Regions of vacuum stability, projected on the heavy scale M vs. electroweak 
scale parameters. All plots are generated at tree level, for fixed value of mh = 115 GeV. 
Left panels: Stability criteria of Eq. (f!5l) . vacuum degeneracy and quantum tunneling are 
denoted by thick gray, thick black and thin labeled lines, respectively. To the left of each 
line the EW vacuum is stable according to that criterion, while to the right it is unstable. 
Right panels: Contours of \fi\ (in GeV) for the same parameter space of the left panels. The 
stability criterion of Eq. ( TT5l) is superimposed to enable easy comparison between the two 
columns. ^ ^ 



a stable vacuum configuration with 120 GeV at tree level. Note that, in order to 

guarantee vacuum stability, it is not a very good practice to assume very heavy uia- The 
reason is that in order for the effective theory to work properly, it is necessary to retain some 
hierarchy between the high scale M and electroweak mass parameters, and so taking any of 
the latter very large can pose a problem. However, positive values for the SUSY breaking 
€2 result in a cancellation in e and so the stability criteria can easily be satisfied through a 
combination of moderate ratios 7714/ ~ 1 and |ei/e| > 1. 

Small quartic coupling limit — \e/e\\ <C 1: An accidental cancellation between the SUSY 
preserving and breaking dimension five operators can lead to |e/ei| <C 1. In this case, for 
negative e±, the RHS of Eq. ( ITBl) approaches zero from above as the ratio |e/ei| decreases, and 
the stability constraint becomes trivially satisfied. A vanishing e corresponds to large €2 > 0.1 
if the LEP bound is satisfied via e\. A finite ratio is then required for the effective expansion 
to remain valid. However, the role of the non-supersymmetric operator €2 is clear. While it 
may not much affect the Higgs mass, this term may be important in stabilizing the potential; 
it thus partially decouples the stability problem from the spectrum. The importance of €2 for 
vacuum stability is illustrated in the lower panels of Figure El We see that turning on even 
a small positive value for €2 rapidly opens up the parameter space corresponding to a stable 
EW vacuum. Note also that the \x parameter is largely unaffected by a change of €2 when 
the Higgs mass is held fixed. This reflects the fact that the stabilization of the potential 
occurs at large field values. 

4.2 Constraints on the scale of new physics 

We now proceed to discuss the constraints imposed by vacuum stability on the scale of new 
physics M. To this end, note that the \i parameter can be traded for M using /1 = e*M/A*. 
We find that imposing stability implies an upper bound on the heavy scale M. For example, 
assuming tra ^ 3mz we obtain from Eq. (fl5l) : 



The upper bound on M for various parameter settings can be read off of Figure [3 Increasing 
ttla, €2 or lowering tan/3 all act to weaken the bound; this behavior is readily understood 
from Eq. ffTTj) . At tan/5 = 5 and with |e| = 0.05, for instance, we obtain |M/Ai| < YItxia- 

Of course, while it is the main issue of the current paper, vacuum stability is not the sole 
source of constraints on the scale of new physics. As an example, in Appendix O we show 
that electroweak precision tests (EWPTs) result in a lower bound for the heavy scale M, 



Put in conjunction with the stability constraints, a lower bound like ( 1181) points towards 
a large value for ttia, sizable SUSY breaking €2, a small tan/3 or some combination of the 
above. Since the dimension six operators responsible for the leading Higgs mass shift and 
for potential stabilization are different from the ones which affect EWPTs, the combined 




(17) 



M > 8 TeV. 



(18) 
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bounds on M may be interpreted in two ways. On the one hand, one may assume generic 
structure for the dimensionless coefficients of all nonrenormalizable BMSSM operators. In 
this case, if quantum corrections are small (e.g., say, both stops are found at the LHC), then 
taking the bound ( Ti~8l) together with the stability constraints leaves a very narrow range for 
the scale of the heavy BMSSM physics. Conversely, one may interpret these two bounds as 
hinting towards some suppression pattern in the microscopic extension behind the BMSSM 
theory. 



4.3 The lightest Higgs mass 

The stability criteria can also be expressed in terms of constraints on the lightest Higgs 
mass. We choose to present the resulting relation between rrih and the heavy scale M. For 
clarity, we consider itla ^ 3m^ and expand to zeroth order in m 2 z /m\. Assuming a negative 
quartic e < 0, we can convert Eq. ( TT7T) into a limit on ei, put in terms of m^tan/?, e 2 and 
M. Plugging the stability constraint in this form into the expression for the Higgs mass, 
Eq. ([7]), we obtain 

m| 2nn 4e 2 sin 2/3(1 + sin 2(3) lQm A sin 2(3^(1 + sin 2/3) /2 

It is illuminating to compare Eq. ( fl9l) with Eq. (171). In particular, notice that in the e 2 term 
a suppression factor of sin 2 2(3 appearing in the spectrum equation gets replaced by a larger 
factor of sin 2(3(1 + sin 2(3) in the stability constraint. This result is explained as follows. 
First, recall that vacuum stability is trivially ensured if the overall quartic coupling along 
the relevant D-flat direction is non- negative, e = e 2 /4 + e\ > 0. Given some value of e 2 , 
substituting t\ = — e 2 /4 (which saturates the relation with a vanishing e) into Eq. (J7J) gives 
the second term on the RHS of ( fl9l) . Hence the last term on the RHS of (fT9l) represents 
the extra gain due to allowing a negative e in the detailed analysis. This gain can indeed 
be significant: for example, at tan (3 = 5(3) the numerical coefficient in front of the ratio 
171a/\M/Xi\ equals ~ 9(16), so that a sizable shift for the Higgs mass is possible even with 
rriA <C M. In addition, if e 2 ~ (in the case where the Beyond-MSSM new physics threshold 
is supersymmetric, for instance) then the last term in Eq. ( fl9l) provides the maximal tree 
level mass shift consistent with vacuum stability. 

The above analysis is illustrated in Figured In the left panel we hold m^,tan/3 and e 2 
fixed. The thick curve depicts the full analytical constraint derived from Eqs. d7J and (fT71) . 
for which Eq. (I19p represents the large limit. (It is easy to verify that Eq. (TIT?]) follows 
this curve to a very good approximation.) In order to keep track of the value of fi, and so of 
ei, we plot contours of constant as thin labeled lines. In the right panel we explore the 
range of accessible in the BMSSM at tree level. We find that with large rriA, small tan (3 
and a sizable e 2 , values of rrih « 140 GeV can be reached at tree level, in a stable vacuum 
configuration, even for |M/Ax| > 8 TeV. Considering the right panel of Figure [3] and using 
Eq. (J7J, we find that |ei| becomes larger than 0.1 at rrih > 148, 130, 113 GeV for the black, 
the dark gray and the light gray curves, respectively. For larger values of rrih the expansion 
to linear order in e\ eventually begins to break down. Thus even if very low values of Mj \\ 
are acceptable, the stability curve cannot be trusted to arbitrarily large m h ; we will touch 
upon this issue again in the next subsection. 
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Figure 3: Left: Stability constraint in the (m^, M) plane, for mA,tan/3 and e 2 as indicated. 
The thick line corresponds to the full analytical approximation (see text). Below this line 
the EW vacuum is stable. Thin contours denote values of Right: Stability constraint 
for different values of 62,rriA and tan/3. 

4.4 Quantum corrections 

The radiative corrections to the scalar potential are dominated by top and stop loops in 
the moderate tan f3 scenario. A (if not The) noteworthy feature of the BMSSM framework 
is that it allows quantum corrections to be small [TT]. As already implied in this work, a 
completely supersymmetric top— stop sector is consistent with both the LEP bound on the 
lightest Higgs mass and stability considerations. Nevertheless, it is interesting to check to 
what extent our results are affected by moderate soft SUSY breaking masses for the stops. 
Here we do not attempt to study the stop sector in full detail. We find it sufficient for our 
purpose to neglect stop mixing and D-terms, adopt a common soft mass for both stops and 
work to one-loop order. Under these simplifications, quantum corrections modify the scalar 
potential by the expression: 



8Vi 



16tt 2 



Q 2 2 J tvr/ v Q 2 2 



(20) 



The field dependent masses are given by m t (4>) = yt<p2, m2 ((p) = m ?(0) + m sto P > with m stop 
the soft stop mass and yt the top Yukawa. Q is the renormalization scale, which we choose to 
be mz- In writing Eq. (1201) . we impose renormalization conditions such that Eq. ([3]) remains 
unchanged. In Figure H] we repeat the stability analysis with quantum corrections, focusing 
for concreteness on the (M, m^) and (M, e 2 ) planes. We learn that soft masses for stops 
stabilize the potential, effectively removing the tree level upper bound on M for mi > 400 
GeV. 

Finally, Figure [5] depicts the effect of quantum corrections on the results derived for the 
lightest Higgs boson mass. Comparing with the right panel of Figure [3] we find, as expected, 
that quantum corrections increase the upper bound on m^. A combination of small BMSSM 
operators e < 0.05 with modest radiative corrections from stops at the 200—250 GeV range 
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can easily and naturally accommodate the experimental bound on m^, even for M > 10 
TeV. Figure includes the tree level numerical stability criterion of the tunneling action, as 
well as the analytical criterion Eq. (jl"9l) . Comparing the two we find that the latter provides 
a reasonable, though slightly conservative approximation to the numerical bound for values 
of mh smaller than ~ 140 GeV. For mh > 140 GeV, |ei| becomes larger than ~ 0.13 and 
the expansion to 0(e) begins to fail, demonstrating the regime of validity of our analytical 
arguments. 



e =0 m =300GeV 

2 A 




MA., [TeV] M/^ [TeV] 



Figure 4: Contours of the tunneling action, including quantum corrections due to top and 
stop loops with a common soft stop mass of 0, 200 and 400 GeV. The allowed parameter 
region lies to the left of the contours. The plots are generated for tan/? = 5 and = 115 
GeV. 
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Figure 5: Contours of the tunneling action, including quantum corrections due to top and 
stop loops with a common soft stop mass of 0, 200 and 400 GeV. The allowed parameter 
region lies below the contours. The plot is generated for fixed m^tan/? and e 2 as indicated. 
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5 Conclusions 



In this paper we analyzed the vacuum structure of the BMSSM Higgs sector. We showed 
that the effective dimension five operators which lift the lightest Higgs mass are potentially 
harmful, as they are capable of destabilizing the scalar potential. It is easy to ensure that 
MSSM D-terms prevent the instability from occurring over most of the field space, for scales 
(f) < M, by imposing e 2 < 1/15. This condition does not exclude a significant shift to the 
Higgs mass. Furthermore, it stands in accordance with the desire to keep the effective theory 
expansion under control. Along the D-flat directions, however, the MSSM D-terms vanish 
and the leading BMSSM correction is dominant. Thus if the quartic coupling along one 
of these directions is negative the effective expansion must be taken beyond leading order. 
Scrutinizing the effective theory to order 1/M 2 , we were able to show that the stability 
of the scalar potential is controlled by only one scaling dimension six operator, which is 
super symmetric and positive definite. This operator is correlated with the supersymmetric 
effective dimension five term which is the most relevant for lifting the lightest Higgs mass. 
Hence the nonrenormalizable part of the theory can cure itself, even though a remote vacuum 
may emerge before the cutoff scale of the effective theory. 

In order to deal with a non-trivial potential structure, a set of criteria was developed from 
which relations between BMSSM parameters were inferred guaranteeing vacuum stability. 
In particular, by analyzing the potential along the D-flat directions we derived an approxi- 
mate analytical criterion, Eq. ffT5"]) . whose robustness was demonstrated by means of a full 
numerical study. Using this criterion we showed, for example, that if the LEP bound on 
the Higgs mass is accommodated at tree level in the BMSSM, then the stability of the EW 
vacuum is ensured provided that htla > Additionally, very low values of tan/3 ~ 0(1) 
are allowed and even favored despite the diminished MSSM contribution to the Higgs mass. 
Interestingly, at the classical level, vacuum stability implies an upper bound on the heavy 
scale M. This bound is better defined than what abstract fine-tuning arguments would 
suggest, and, since the experimental data can be accommodated with small fi/M < 0.05, 
is also typically stronger. Put in conjunction with generic lower bounds on M, arising for 
instance from electroweak precision tests, the analysis either highly constrains the BMSSM 
parameter space or directs us towards a non-generic coupling structure in the effective theory. 

Vacuum stability also dictates an upper bound on the lightest Higgs mass accessible via 
the leading effective operators. However, for M in the few TeV range, > 140 GeV can 
still be accommodated at tree level, as a result of the potential stabilization provided by 
the supersymmetric dimension six operator. While a completely supersymmetric top— stop 
sector is allowed in the BMSSM (as far as the experimental bound on the lightest Higgs 
mass is regarded), quantum corrections due to stop loops are effective in stabilizing the 
potential, and two stop states at ~ 300 GeV suffice to significantly weaken the constraints 
due to vacuum stability. Indeed, light stops just around the corner for upcoming colliders 
can provide a relatively unimportant direct contribution for the Higgs spectrum, yet be 
significant for vacuum stabilization at large field values. 

Finding both stops not too far above m t would already be a smoking gun for a Beyond- 
MSSM extension. Information on the spectrum of charginos, neutralinos and the Higgs 
scalars, together with stability considerations, could then be used in order to extract addi- 
tional constraints on dimension five operators and low energy MSSM parameters, as well as 
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put an upper bound on the heavy scale M. 
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A CP violating and charge breaking background fields 

Under the assumption of real parameters, the full effective potential including \ and p fields 
is given by 



+2e 2 



2 2 
^2 „,2\ j£ i 9z | ,,2 i „2 j,2\ 2 i 9 A 2, ,2 



x 2 )<t>l + f(rt + x 2 + P 2 -<t>lY + y Y 



+4 



02 (01 + X 2 + p 2 + 02) (0 2 + X 2 )- 



First derivatives of interest follow 
dV 



dx 2 



ml - 4 ei 0!0 2 - 2e 2 2 + if. U\ + x * + p 2 



+4 



^ (20? + 2 X 2 + ^ + p 2) 



(21) 



(22) 



91/ 



4 



D+4 



>2M + X 2 ) + y<& (23) 



Since only the squares of x an d p appear in the potential, a trivial extremum solution with 
X = p = always exists. CP violating (CPV) and charge breaking (CB) extrema require 
nontrivial solutions for a vanishing RHS in Eqs. ( 1221) and ( 1231) . respectively. 

We first derive the condition ensuring that the EW vacuum is CP and charge conserving. 
Demanding that the potential has a minimum at <f>i — v cos ft, 2 — v sin ft and p = x — 
and expanding to leading order in e, we obtain that an additional, nontrivial solution with 
0i = v cos ft, 2 = f sin ft and x 7^ is never allowed, while a solution with p ^ requires 
p 2 ^ _4 s i n 2 ^ {m 2 A + m 2 ^ + 2e 2 f 2 ) / g 2 z . Hence x — P — is the only consistent solution at 
the EW vacuum, as long as: 



m A + m w + 2e 2V 2 > 0. 



(24) 



Considering Eq. ([9]) and the validity of the effective theory expansion, we find that this 
relation is always satisfied for itla > mz- 

We now move on to find out how nonzero x or /and p affect the potential at large field 
values. In particular, we are interested in configurations of x, P which extremize the potential 
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in the vicinity of the tunneling path. Finding such extrema is important, since they may 
shift the saddle point and/or the remote vacuum from the (0i,02) plane and so alter the 
numerical computation of the stability criteria. Note that for large fields, unlike in the 
local neighborhood of the EW vacuum, one must keep the scaling dimension six operator 
appearing in the potential ( |2Ti) . In addition to the trivial solution {x = 0,p = 0}, there are 
three exclusive possibilities: {x / 0,p = 0}, {x = 0, p 7^ 0} or {x 7^ 0, p 7^ 0}. To find 
out which of them is consistent, we rewrite the condition of vanishing right hand sides in 
Eqs. (J22D and (J2HD as 

(25) 



X 2 = A C pv(<Pi, 02) + B C p V {<t>\, 02)p 2 , P 2 = A CB {(j>u 02) + B CB {(j>i, h)x 2 



Solutions with {x 7^ 0,p = 0} or {x = 0,p 7^ 0} only exist for Acpv > or Aqb > 
respectively. Moreover, since the second derivatives of the potential with respect to x 2 and 
p 2 are positive definite, whenever one of these solutions becomes available within some region 
of the (0i, 02) plane it will always minimize the potential energy, rendering the trivial solution 
unfavorable. Hence if the tunneling path traverses these regions, CPV or CB will turn on. A 
solution with {x 7^ 0,p 7^ 0} requires that both x 2 — (Acpv — AcbBcpv) I '(1 — BcbBcpv) 
and p 2 = (Aqb — AcpvBcb)/ {1 — BcbBcpv) be positive. If such a solution exists, it may 
or may not become energetically favored over solutions with only x 01 P nonzero. 

We now describe how the x and p fields are dealt with in the numerical procedures 
presented in Section 13.21 The criterion of vacuum degeneracy requires computing the 
value of the potential at the remote minimum. For the tunneling action, according to the 
approximation we adopt (see Appendix [B]) , the location of and potential value at both the 
remote minimum and saddle point are needed. To find these extrema in the presence of 
possible CPV and CB configurations, we numerically evaluate the full potential over the 
(0i,02) plane where x and p are replaced by the four classes of solutions defined above. 
Then we retain all consistent solutions extremizing the potential along both the 0i and 02 
directions. While the remote minimum is unique, as it has to minimize the potential along 
all directions, there may be several saddle points corresponding to the different classes of 
solutions for x and p. We compute the tunneling action assuming the path passes through 
the saddle point of minimum potential energy. 

After the dust settles it turns out that in typical cases in which the EW vacuum is found 
to be stable, no CPV or CB arise in regions of field space that are relevant for the analysis, 
or otherwise the effect is very small. Thus the discussion conducted in Section I3TT1 where for 
clarity we have set x = P = 0, holds true up to minor changes. In order to see how this result 
comes about, it is useful once again to consider the problem close to the D-flat directions. 
Since away from them positive D-terms raise the potential, the tunneling path is confined 
to remain near this region. Thus removing CPV and CB from the flat directions protects 
the saddle point and the remote vacuum from developing nonzero x or P- We can derive 
simple analytical criteria to exclude the formation of CPV or CB extrema along the D-flat 
directions. To this end, let us consider e x < in which case potential instability occurs for 
negative 0i. Setting — 0i = 2 = <fi/V2 and e = €2/4 + e\ < 0, the extremum equations 
reduce to: 



9l 



+ 4 




2 p 2 , (26) 
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We find that CPV and CB are banned along the D-flat directions provided the following 
conditions hold, respectively: 



,2 



up to 0(m 2 z /m 2 A ) corrections, where /(/3) = (1 + sin 2j3)/ sin 2 /? satisfies 1 < /(/?) < 4 
for tan/5 > 1. These two conditions can be compared to the approximate analytical 
criterion derived in Section [37X1 which ensures that the EW vacuum is stable, namely m\(l + 
sin 2/5) / 2 > 2(e/ei) 2 . Regarding CB, we find that vacuum stability guarantees no CB 
turns on along the D-flat direction provided that |ei| < 0.2. In view of Figured! as well as 
considering Eq. @ and the validity of the effective theory expansion, it is clear that this 
relation does not constitute a real compromise for the size of e\. The appearance of CPV is 
a less marginal effect. However, the form of the inequalities ( |28|) and the stability constraint 
implies that the stable parameter space exhibits little CPV; numerically we find that in most 
stable scenarios CPV does not occur along the tunneling path, or else does not encompass 
the saddle point nor the remote minimum. We conclude that imposing vacuum stability 
typically renders CB and CPV configurations irrelevant in the calculation. Note that since 
the tunneling path does not pass strictly along, but merely in the vicinity of the D-flat 
direction, the above analysis provides an intuitive argument only. However, as previously 
mentioned, numerically scanning the relevant parameter space we find that in practice it is 
indeed safe to neglect \ and p in the analysis of vacuum stability in the BMSSM with real 
Lagrangian parameters. 



B Quantum tunneling 

We are interested in analyzing the stability of configurations wherein the potential exhibits 
a remote vacuum, in the presence of which our EW vacuum may be metastable. At zero 
temperature, one should compute the rate of quantum tunneling from the EW to the remote 
vacuum. For completeness, let us briefly review the theoretical set up before going into the 
details of our implementation. In a semi-classical approach the tunneling rate (per unit time 
and volume) is given by the WKB approximation [22] : 

r ~ Aexp (-B) , B = S E [<p b ] - ^[0 false ], (29) 

where Se[4>] is the Euclidean action and A is the determinant of the Gaussian fluctuations 
around 0b, the so-called bounce solution of the equation of motion. The precise value of the 
prefactor A plays a minor role in comparison to the exponential suppression factor. Hence 
we make the conservative assumption that the prefactor saturates the highest scale in the 
BMSSM framework: A ~ M 4 . The bounce solution is 0(4) symmetric [23], i.e. is only a 
function of r = \J x 2 + t 2 E . It solves the equation of motion 

d 2 (f) b 3 d(f) b 

~dr 2 r~dr ~ V ^ = ' ( 30 ) 
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subject to the "bouncing" boundary conditions: (j)' b (0) = and 0b(oo) = 0f a i se . 

In configuration space, quantum tunneling proceeds by nucleation, growth and percola- 
tion of true vacuum bubbles surrounded by a metastable environment. At a given time t, 
the portion V(t) of the volume of the universe filled by true vacuum bubbles is controlled 
by the tunneling rate and the expansion of the universe through the following relation [21] 



V(t) = 1 - exp 



Air C % 

-| ^ dt'mdu?) , (si) 



where dn{t) ~ t is the horizon distance. For the case of quantum tunneling (T = 0), the rate 
T, given by ([2"§]) . is time independent. Hence today, for t ~ 10 10 yr ~ e 100 ?; -1 , the portion 
of the universe in the stable phase is given approximately by: 



V(t ) ~ 1 — exp 



-O(lO)cxy) ( 41og^ - B j (>/,,,; 



(32) 



Therefore to ensure that the metastable vacuum has a life time longer than the present age of 
the universe, or equivalently that V (to) — 0, it is enough to impose the following constraint 
on the bounce action: 

5 > 400 + 4 log— , (33) 

v 

where the second term, coming from the exponential prefactor A, contributes at most 0(10%) 
for M in the TeV range. 

In a one-dimensional field space (as in the Standard Model), the correct bounce solution 
is easily found by use of a simple "shooting" numerical method. However, in the case of 
multiple fields, the search of the bounce solution is a non-trivial numerical task. We make 
use of several simplifying assumptions, justified in the BMSSM setup under consideration. 

The first simplification consists of projecting the action on a given path in field space, 
which is known to be close to the correct tunneling trajectory. This path is fully characterized 
by a single curvilinear coordinate <p s , i.e. <fi = 0(0 S )- The initial task is then greatly 
simplified as the problem has been reduced to a one dimensional search for the bounce 
solution. However such a projection is in general hard to estimate. One exception is when the 
potential extrema are close to alignment in field space. In this case the tunneling trajectory 
may be well approximated by a straight line joining the true and false vacua [25J. Projecting 
on such a path allows one to easily compute the bounce action using a traditional shooting 
method. We stress that this approach is consistent as the approximate alignment that it 
requires typically arises in the BMSSM. Indeed, as is apparent in Figure EJ both the remote 
vacuum and the saddle point separating it from the EW vacuum form close to the D-flat 
direction. 

Second, we adopt a triangular approximation for the potential barrier, allowing us to use 
the clear prescription given in [26J. Thus we reserve the more cumbersome shooting method 
only as a check on the triangle calculation, verifying that the two methods agree at the level 
of (9(10%). In Figure [6] we demonstrate the use of the triangle approximation in computing 
the tunneling action. 
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Figure 6: Left: Contour plot of the potential for p = tua = 300 GeV, tan (3 = 2, €2 = 
and ei ~ —0.05, corresponding to tree level Higgs mass of rrih ~ 115 GeV. The straight line 
connects the two minima. Right: Potential projected on the straight line joining the two 
vacua (solid) and the corresponding triangle approximation (dashed) used to estimate the 
tunneling action. <ft s is the coordinate along the straight line. On both plots <fi and V e s are 
in units of v and v 4 respectively. 

C Electroweak constraints 

As we have seen, stability of the scalar potential, arising from effective dimension six 
operators, gives an upper bound on the scale of new physics. Other effective dimension 
six operators involving the Higgs sector also affect the gauge terms in the Lagrangian, 
introducing mass shifts in the gauge sector as well as kinetic mixings. They are therefore 
constrained by electroweak data, in particular the precision electroweak variables S and T 

S = -0.10 ±0.10 , T = -0.08 ±0.11 (34) 

for rrih ~ 115 GeV. There are several effective dimension six operators involving the Higgs 
sector that can contribute to electroweak observables. For example, operators of the form 
aij(Hje v Hi)(Hje v Hj) /M 2 with i,j = u, d in the Kahler potential contribute to the deviation 
of the p parameter from unity by [10] 



/ \ 2 

P - 1 = -4 ( j [a uu sin 4 (3 + a dd cos 4 (3 - a ud cos 2 (3 sin 2 (3] ~ aT. (35) 

For generic coefficients ciij ~ 0(1), the expression in square brackets in ( 135]) is not particularly 
suppressed for any specific value of tan (3. We obtain a lower bound on the heavy scale, 

M>7.7TeV. (36) 
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Regarding the S parameter, consider for instance a superpotential operator of the form 
a WB W a B a H u H d /M 2 , where W a {B a ) denotes the SU(2) L (U(1) Y ) gauge superfield strength. 
Such an operator generates a contribution to S of 

327r t; 2 sin2/j 

S= ^^2M^ aWB - (37) 



Taking awB ~ we obtain a comparable bound, M > 8vsm2^ TeV. 
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